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SUMMAEY 

As a result of a fundamental investigation of the meteorological 
conditions conducive to the formation of Ice on aircraft and a study 
of “the~proce S 8 “of“ alrfoli" thermal— ice“preve nt ion,— previousl-y-der ived— 
equations for calculating the rate of heat transfer from airfoils in 
icing conditions were verified. Knowledge of the manner In which 
water is deposited on and evaporated from the surface of a heated 
airfoil was expanded sufficiently to allow reasonably accurate calcu^ 
lations of airfoil' heat requirements. The research consisted of 
fll^t tests in natural-icing conditions with two 8 -foot-chord heated 
airfoils of different sections. Measurements of the meteorological 
variables conducive to ice formation were made simultaneotisly with the 
procxurement of airfoil thermal data. 

It was concluded that the extent of knowledge on the meteorology 
of icing, the ir^jlngement of water drops on airfoil surfaces, and the 
processes of heat transfer and evaporation from a wetted airfoil 
surface has been increased to a point where the design of heated wings 
on a fundamental, wet-air basis now can be undertaken with reasonable* 
certainty. 


INTRCDUCTION 

For a period of several years, the National Advisory Committee 
for Aeronautics has conducted reseeu*ch on the prevention of ice 
formation on aircraft through the use of heat. During this time, 
reseeu'ch of a fundamental nature on the problem of thermal ice 
prevention was retarded by the more urgent need for development of 
ice-prevention systems for speolflc airplanes in military service. 
Satisfactory wing- and tall-evjrface thermal, ice— prevention systems 
for a Lockheed 12-A, Consolidated B-2k, Boeing B-I 7 , and 
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Curtiss-Wright 0-U6 airplanes (references 1, 2, 3, and U, respectively) 
were designed, fabricated, and tested in natural-icing conditions. 
Windshield thermal ice— prevention systems idiich proved adequate in 
the icing conditions encountered were provided for the 12-A, B-24, 
and C-46 airplanes. Each wing- and tail-surface design was based on 
establishing, for fli^t in clear-air conditions, a surface- 
temperattoa rise above free-stream tenqjerature i^ich experience in 
simulated— and natural— icing conditions had shown to be adequate for 
ice prevention. This eiiq)irical method, while proving satisfactory 
for the airplanes tested, was limited, since it was not established 
on a fundamental basis, and a more baaic procedure foioided on 
designing for the conditions existing in icing clouds was needed. 

Ihe NACA at present is engaged in an investigation to provide a 
fu ndame ntal imderstandlng of the process of thermal ice prevention 
in order (1) to establish a basis for the extrapolation of present 
limited data on heat requirements to meteorological and fll^t condi- 
tions for which test data are not available, (2) to provide data for 
lBq)rovlng the efficiency of thermal ice-prevention equipment, and 
(3) to provide a wet-air, or meteorological, basis for the preparation 
of design specifications for thermal ice— prevention equipment. The 
research consists of an investigation of the zneteorologlcal factors 
conducive to icing, and a study of the heat-transfer "processes which 
govern the operation of thermal ice— prevention equipnmnt for airfoils 
and for windshield configurations. 


The airfoil heat— transfer phase of this investigation consisted 
of the measurement of the factors affecting the transfer. of heat 
from airfoil surfaces during flight in natural-icing conditions. 
These data are correlated with the simultaneous measurements of the 
meteorological parameters which influence the heat-transfer process, 
and are ai3alyzed for the purpose of establishing a wet-air ice- 
prevention design basis for airfoils. 

The first approach to the icing heat— tremsfer problem on a 
fu n da m entetl basis was made in England by Hardy and Mann prior to 
19^ • In this study, a method for the calculation of heat transfer 
from a heated surface subjected to icing conditions was presen'tod 
and substantiated by meastirements in an icing tunnel. Later work by 
Hardy in \dilch these heat-transfer equations were modified for 
general application is presented in references 5 and 6. Reference 5 
contains information on the p>rotectlon of a- 1 i aircraft conponents 
agalxist ice accretion. Reference 6, prepared during a period of 
active participation by tfr. Hardy in the MCA icing research program, 
presents an analysis of the dissipation of heat in conditions of 
icing from a section of the heated wing of the C--h6 airplane 
(reference 4) . 
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Other research In the present MCA inrestigatlon has been 
reported In references 7> 8# 9* Reference 7 gives the first 

measurements in this program of the liguld-^rater concentration in 
clouds. References 8 and 9 deal with the meteorological aspects Of 
icing conditions in strat^ clouds and in precipitation areas of the 
warm-front type. 

Research on the problem of heat transfer from airfoils in condi- 
tions of icing has also been conducted by other laboratories. In 
reference 10, the transfer of heat from surfaces subjected to icing 
conditions on Mount Washington has been studied. The General Electric 
Research Laboratory has conducted a number of investigations on this 
phase of icing. A summary of this work emd a list of imports is 
presented in reference 11. A con 5 )rehenelve report by the Army Air 
Forces on the developmant and application of heated vlngs is contained 
in reference 12. 

In continuation of the present icing ~progr^,~tKe“C^6^IrTli^~^' 
was equipped with special meteorological and elecWlcally heated test 
apparatus, and flown in natinal— icing conditions dixrihg the winters 
of 19k^-ko and 19^6-47. Fll^t tests were conducted mainly along 
airline routes over most of the Ihiited States. The meteorological 
data recorded during the icing conditions encountered in the two 
seasons are presented and disctxssed in references I 3 and l4. 

This report presents an analysis of the data obtained d\n*lng the 
1945—46 and 1946-47 winter seasons with two electrically heated air- 
foil sections. The data were analyzed lusing the heat-transfer 
equations developed by Hardy. (See references 5 and 6.) A considera<- 
tion of the area and rate of water impingement on one of the airfoil 
sections based on an analytical study of water-drop trajectories 
(reference I 5 ) is also presented. An attempt is made to further the 
knowledge of the process .of airfoil thermal ice— prevention. 

The appreciation of the MCA is extended to IMlted Air Lines, 

Inc., the Ikiited States Weather Bui^au, and to the Air Materiel 
Command of the Army Air Forces for aid and cooperation in the research. 
In particular, the services of Major James L. ^ftIrray of the Air 
Materiel Command, Army Air Forces, and Captain Carl M, Christenson 
and First Officer Lyle W, Reynolds of United Air Lines, Mho served as 
pilots of the research airplane, were a valuable aid to the conduct of 
the investigation. 
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!Qie foUovlng nGmenclature I0 uirod throiighout Wls report: 
a radius of water drop, feet 


c 

Cp 


°Pw 


C 

e 


E 

8 

h 

H 

J 

k 

K 

Lb 


airfoil chard length, feet 

specific heat of air at constant pressure, Btu per pound, 
degree Fahrenheit 

specific heat of water at constant pressixre, equal to 1 Btu 
per poimd, degree Fahrenheit 

concentration factor, defined in equation (6), dimensionless 

saturation vapor pressure with respect to water, millimeters 
of mercury 

wate3>-drop collection efficiency, defined in equation (lO) 

acceleration of gravity, equal to 32*2 feet per second, second 

convective surface heat-transfer coefficient, Btu per hour, 
square foot, degree Fahrenheit 

total sxnrface heat-transfer coefficient, Btu per ho;n:, square 
foot, degree Fahreideit 

mechanical equiveJ-ent of heat, equal to 778 foot-pounds per 
Btu 

thermal conductivity, Btu per second, square foot, degree 
Fahrenheit per foot 

dimensionless drop-inertia quantity, defined in equation (^) 

latent heat of vaporization of water at stu'face ten^ierature, 
Btu per pound 

liquid-water concentration of icing cloud, poimds of water 
per cubic foot of air 


1% wei^t rate of water-drop liapingement per unit of surface 
eurea, poimds per ho^, square foot 

Mb weight rate of water flow aft of area of \mter-drop iBq)lnge- 
ment per foot of span for one side of airfoil, pounds per 
hour, foot 
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weight rate of water-drop in^jingement per foot of span for one 
side of airfoil, poxmds per hover, foot 

concentration of liquid water contained in drops of each size 
in a drop-size distrihution, pounds of water per cubic foot 
of air 

barometric pressure, millimeters of mercury 

Prandtl number (cp|a/k), dimensionless 

tmit rate of heat flow, Btu per hour, square foot 

Eeynolds number for airfoil (Vey/ti), dimensionless 

f ree— st ream Reynolds number of water drop relative to speed of 
airfoil ^(SVaT/tr)'; dlmeiieionless 

distance measured chordwise along airfoil surface from staigna- 
tion point, feet 

tenq)erature, degrees Fahrenheit 

local velocity Just outside boundary layer, feet per second 
free— stream velocity, feet per second 

weight rate of evaporation of water per unit of surface area, 
pounds per hour, square foot 

weight rate of evaporation of water per foot of span for one 
side of airfoil, pounds per hour, foot 

distance measvired chordwise along airfoil chord line from 
zero— percent chord point, feet 

evaporation factor, defined in equation (22), dimensionless 

airfoil ordinate, feet 

starting distance of water drop above projected chord line of 
airfoil, feet 

pressvire altitude, feet 

ratio of saturated to dry adiabatic lapse rates 

exponent of Prandtl munber, 1/2 for laminar flow, I /3 for 
turbulent flow / 
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7 epeclfic wel^t of air, pounds per cubic foot 

7v specific wel^t of water, eg.ua! to 62,4 pounds per cubic foot 

|i viscosity of air, pounds per second, foot 

Subscripts 

1 refers to conditions at ed^e of boundary layer 

k kinetic' 

o refers to free-stream conditions 

me mean effective 

s refers to conditions at airfoil surface 

SL sea level 


ANALYSIS 

During flight in icing conditions a heated wing is cooled by 
convective heat transfer, by evaposration of the water on Jie surface, 
and, in the region of droplet interception, by the water striking 
the wing.i The rate at which heat must be supplied in order to 
maintain the wing surface at a specified temperature is, therefore, a 
function of the rates of convection, evaporation, water ir^jingement. 
Equations for expressing this heat requirement are presented in 
references 5 and 6, These equations, with sli^t modification, are 
used throu^out this report. 

Expressed as an equation, the unit heat loss q from a partially 
or conq)letely wetted surface exposed to icing conditions may be stated: 

<1 = 4w + Ic + <le ( 1 ) 

where 

qw heat loss due to warming the Intercepted water 

qc heat loss due to forced convection 

q^ heat loss due to evaporation of the Implnglisg water 

Each of these individual heat flows will be analyzed. 


^The heat loss due to radiation is small and can be neglected< 
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Heat Loss Dus to Wamlog the Intercepted Water 

In the region vtnsre vater droplets strike the wing, the heat 
required per unit area to heat the water to surface temperature Is 

4w “ “ ^'*'0 j 

The term Atj^^ Is the kinetic tenqperatiire rise of the water caused 
hy stoppage of the droplets as they strike the wing. The value of 
Atfci, Is given hy 


Atky 


yg 

SgJCpy. 


( 3 ) 


where V Is^tfi^free-s't^aitt vel^lty In feet per second. The value 
of Atky Is less than 2® Fahx*enhelt for airplane speeds up to 200 
miles an hour and, for the calcxilatlons presented In this report, the 
term has heen neglected. Equation (2) thus 'becomes: 


4w ® Ha (■fcer-'fco) 


( 4 ) 


Ihe wel£^t rate of water Impingement on the wing, the area of 
impingement, and the distrl'butlon of the water over that area are 
important factors In the heat— transfer analysis. In addition to the 
effect of the amount of water intercepted on the value of q^ in 
equation (U), the evaluation of % provides an indication of the 
quantity of water idiich mmst "be maintained in a liqviid state until 
it either evaporates or runs off the trailing edge if the formation 
of ice aft of the area of impingement, normally termed "runback," is 
to he avoided. Hie area of impingement influences the extent of 
heated region to he provided at the leading edge, while knowledge of 
the distribution of water impingement is required in the calculation 
of the heating requirement in areas where water is striking. 

Calculations have heen made hy Glauert (reference I6) for the 
trajectories of water drops about cylinders and an airfoil. In this 
work the assumption vas made that the drops obeyed Stokes* law of 
resistance. At the speeds of flight, however, Stokes* law no longer 
strictly holds, and Langmmir and Blodgett (reference I7) computed a 
series of drop trajectories about cylinders, spheres, and ribbons, 
taking into consideration deviations from Stokes* law. These compu- 
tations were undertaken on the assumption that the trajectories for 
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cylinders would apply to airfoils if the airfoil were replaced by an 
"equivalent" cylinder (reference 12), ; 

Preliminary calculations based on references l6 and I 7 indicated 
that, for large values of drop size and airspeed, the assmq)tion of 
the equivalent cylinder would not hold for airfoils. Therefore, more 
extensive calculations were undertaken to determine the drop trajec- 
tories for one of the test airfoils of this research, an NACA 0012 
airfoil at 0 angle of attack. In these calculations, presented in 
detail in reference 15, a Joukowski airfoil (the contoiir of which 
closely approximates that of the KACA 0012) was used to supply the. 
stream lines since the Joukowski stream lines and velocity field can 
be computed with relative ease. The basic equations presented in 
•reference I 6 were used with modifications for deviation from Stokes' 
law as given in reference 17» The procedure followed was to start a 
given distance forward of the airfoil and calculate the paths of the 
drops using a step-by— step integration process, Eesults of these 
computations are presented in figure 1, Ihe cxjrves shown establish 
the distance s, measured from the stagnation point, at which a given 
drop will strike the airfoil when starting a distance y© above 
the projected chord line. Ciirves are presented for various values 
of K, where 


It should be noted that the curves of figure 1 apply strictly only 
for a drop Reynolds number Eu of 95. that is, only for partic- 
ular combinations of drop size, airspeed, altitude, and air ten^jera— 
tvire. The value of 95.^5 was chosen as being the Reynolds nmber 
corresponding to average conditions of drop size, airspeed, altitude, 
and air tenq>eratvjre experienced during the tests of this investiga- 
tion. However, the curves of figure 1 can be used for a range of 
Reynolds numbers on either side of 95.65 without serious error. Due 
to practical consid.eratioixs, these curves were used in the analysis 
of the data presented in this report, even though the Reynolds 
number differed somewhat for every case. 

Area of water impingement,— The end points of the curves shown 
in figure 1 denote the extreme location at which drops of a partic- 
ular K value will strike the airfoil. Beyond this value of s/c no 
drops of this K value will hit. Thus, the broken line in figure 1 
establishes the area of impingement for all values of K, 

Rate of water impingement .— The rate of water inq>ingement at a 
specified point on an airfoil is a function of the area of impingement, 
the velocity of flight, the liquid— water concentration of the air 
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stream^ and the dietrlhution of the intercepted vater mer the surface* 
»Phi H latter factor, called the concentration factor C is represented 
l>7 the ratio of to s, or: 


For point Yalnes, 


C 


Ayo 


or more exactly, 

0 - ^ ( 6 ) 

ds 

The welj^t rate of-wter iiiq^iiigement per uxilt of^surface area-1^ 

pouxids per hotcr> square foot^ then^ Is 


» 3600 Y m C (7) 

It is apparent from equation (6) that C is simply the slope of 
the curves shom in figure 1* A plot of the measured slopes of these 
curves as a function of s/c is presented in figure 2* Using values 
of C ohtalz»d from figure 2, the 'uel£^t rate of vater liiq>lngement 
at any point on the surface can he calcxilated from equation (7). 

In the case of a cloud, vhex« the vater drops are not of uniform 
size, hut Instead foUov a pattezn of size dlstrlhution, the rate of 
is^lngement can he canq>uted if the dlstrlhution is hnovn or assumed* 
The rate of vater impingement at any point is the" sum of all the rates 
of Is^ingement of the Tolume of vater contained in each drop size* 
Equation (T) then hecooDoes 

nC (8) 

where n is the concentration of liquid vater contcdned in drops of 
a particular size and C is the concentration factor for the K 
value corresponding to that drop size* ,, 

In order to estahllsh the posslhlllty of runhack f otmlng aft of 
the heated area of a vlng, it is necessary to know the total quantity 
of vater Intercepted per unit of vlng span* This rate of impingement, 
denoted as He In pounds per hour, foot span, is given hy 


Ma » 3600 Y 


I 
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Mg 


/' 


Ifg, dS 


(9) 


A more rapid method for the eraluation of Mg utilizes a curve of 
collection efficiency E as a function of K (fig. 3)« Collection 
efficiency is defined as 


j. ^ ypumit 

ymax 


( 10 ) 


where yoiimit value of y© for >diich drops of a particular 

K value Just miss the airfoil, and ymax is the maximum ordinate of 
the airfoil. The equation for computing Mg, then, is 


» 3600 EVm ymax (H) 

Using figure 3, the rate of water impingement can he ccoputed for 
each of the drop sizes in the assumed or measured drop-size distri- 
bution. The total rate of iaq)ingement is the summation of these 
individual rates. 


Heat Loss Due to Forced Comrecticui 

3he unit heat flow from the surface of a body in an air stream 
resulting from convective heat transfer can be expressed: 


qc = h (tg— toic) 


( 12 ) 


where tg is the surface tenq)erature and tofc is the kinetic tenqpei'- 
ature of the fTee-stream air at the point for which the heat flow is 
being conqjuted. The factor h is the. convective heat-transfer 
cpefficient and may be evaluated by measurements in clear air or by 
calculation using the methods presented in references I8 and I9. 
Evaluation of the term tgj^ will now be discussed. ^ 

The surface of an unheated wing moving through the air will 
assume a temperature somewhat hi^er than that of the free air stream 
because of stoppage of the air particles in the boundary layer next 
to the surface. ■ ^is temperature rise is of importance in the calcic 
lation of heat requirements for ice prevention in that it establishes 
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the datum point from which the temperature of the surface must he 
raised, to obtain the desired tenperature, tg.. The value of the 
teu^erature rise in clear air, from equations derived in reference 5, 
is, for laminar flow. 


Atjj. = -2^ r 1 
^ 2gJcp L 


and for turbulent flow. 


Fl 

ggJc, L 


where U is the local velocity Just outside the boundary layer at the 
point along the surface where the value of At]^ is being calculated. 

In clouds, the kinetic temperat\u:*e rise is reduced, due to 
evaporation of water fi'om the surface. Assuming the surface is 
conq>letely wetted with water, the value of the tenq^rature rise for 
laminar flow becomes 




where 


01 “ e© ^ (16) 

* o 


(13) 

(I**) 


and eot Is the vapor pressure at saturation at the wet kinetic tem- 
perature, toif. Fn»e value of toj; is 


or 


tojj “ tg + Atj. 
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Sqviatlon (l?) Is for lamina-y flow, Hie eq.iaation for tiir'bulent flow 
Is the same, hut with the exponent of ?r changed to I/ 3 . It can he 
seen that this equation must he solved hj trleil, since the value of 
eo^ Is dependent upon the tenqeratiire to]^;, 

laqperlments In clouds. In the process of calibrating a free-alr 
thermometer Installation (reference I 3 ), showed that hy multiplying 
the cleaj>^lr kinetic— terq^rature rise hy the ratio of the saturated 
to the dry adlahatlc lapse rates, good agreement between the values 
of kinetic temperature rise calculated In this manner and the measured 
values was obtained. Since use of the ratio of the adlahatlc lapse 
rates was substantiated experimentally, and since equation (I 7 ) must 
he solved hy trial, a somewhat laborious procedure, the foUo^ng 
equations were iised in this report to calculate values of to]^: 

For laminar flow. 


^°k “ 2iJcJ 


and for turbulent flow. 

Values of as/od> the ratio of the wet- to the dry-adlahatlc lapse 
rates, are obtained from figure 4, Hie use of the lapse-rate ratio 
in equations (I 8 ) and (I 9 ) is seml-enq)lrical. Hie limitations of 
this slsq)llflcation In the calculation of kinetic— teoqieratuz^ rise of 
airfoil surfaces In clouds are not known. Below speeds of 200 miles 
per hour, however, these equations can he used with small error, since 
the kinetic— temperature rise Is low. 


( 18 ) 

[l-p (19) 


Beat Loss Due to Evaporation of the V&ter on the Surface 

Hie amount of heat removed from a wetted surface as a result of 
the evaporation of water on that surface can he expressed: 
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iVom. reference 6 the relation hetveen qe and the convective 
heat^transfer coefficient h can he expressed for a con^etely 
vetted surface as: 


Oe =h (X -1) (tg-toj^) . (21) 



By substituting average values for Ls and cp^ equation (22) can 
he revrltten 




Xhe values chos^ for Ls and cp are 1100 Btu per pound and 0.24 
Btu per pound, 'Tahrenhelt, respectively* Ohe factor PslAi is the 
ratio of the standard sear-level prese:n‘e to the local static pressure. 

It should he noted that the evaporation factor X applies only 
when the surface Is cos^letely vetted. If only psirtlal vetness 
prevails, the value of X must he modified according to the degree 
of vetness* 


Total Heat Loss from a Wetted Svirface 

Summarizing the heat losses due to vater lIl^>lng^m9nt, convection, 
and evaporation, equation (l) can he vrltten: 

\ 

q. = Ha (t|^to) + h (ts— to]f) + h (X — l) (tg— tQjj.) 


vhlch reduces to 


(24) 


q, = Ma (i'S—'to) + h X (ts— toj^) 
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Aft of the region of vater iii^>inge!DBnt, as 0 and equation (2U) 
hecomes 


q = h X (ts-to^) (25) 


DEScEipnaBr of equipmemt 

All tests reported herein were made In the 0-46 airplane shown 
In figure 5. The airplane had teen modified to provide thermal Ice- 
preventlon equipment for wings, enqpennage, windshield, and propellers. 
A description of the thermal system for the wings wnfl empennage Is 
given In reference 20. The windshield system was altered for the 
flights as described In reference 21. Protection for the propellers 
was provided by electrically heated blade shoes. 

I 

The meteorological equipment used during the tests to measure 
the free-air temperature, liquid-water concentration, drop size and 
drop-size distribution Is described In references 13 l4. 

Two electrically heated test airfoils were used to obtain 
fundamental data on the process of wing thermal Ice prevention. Each 
airfoil was mounted vertically on top of the fuselage of the C-46 
airplane, as shown In figure 5. The tost airfoil Installed during 
the winter of 1945-46 had an HACA 0012 section. For the tests in the 
winter of 1946-47 the airfoil had an HACA 65,2-016 section in order 
to provide test data for low-drag sections, as well as conventional 
sections. Both sections are symmetrical, and the models were 
Installed with the chord line in the plane of symmetry of the alj>- 
plane; that is, at zero aiigle of attack for unyawed fllc^t. Ordinates 
for an HACA 0012 airfoil are. given In reference 22, and for an 
HACA 65,2— 0l6 airfoil. In reference 23. Figure 6(a) shows the 
HACA 0012 airfoil mounted on the fuselage, the HACA 65,2-016 
airfoil was mounted as shown in figiure 6(b). A clear plastic 
blister, shown in figures 6(a) and 6(b), allowed the airfoils to 
be viewed and photographed In fll^t. 

Both airfoils had an 8-foot chord and a 4.7-foot span, with a 
faired sqxiare tip. A heated test section of 1-foot speua was located 
2 feet above the top of the fuselage. It had been determined 
previously, by means of a pressure survey, tliat the test-section 
location was well above the edge of the fuselage boundary layer. 
Electrically heated guard sections were built around the leading- 
edge region on both sides of the test section for the purpose of 
preventing any dietxa*bance of the air flowing over the test section 
which mi^t have been caused by Ice accretions in the region of the 
guard sections. 
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NACA 65,2-016 Airfoil Model 

Conatrtictlon details ,of the 65,2-016 airfoil model are shown In 
figure 7. The metal portlcm of the structio'e consisted of alumlnwn 
rlhs and skin supported from the fuselage by two spars. The test 
section was made up of a 3/8-lnchr-thlck plastic base and a sheet of 
plastic— Impregnated fabric, l/61»— Inch-thick, on top of which 1/2— Inch 
wide, 0.002— Inch-thick, electrical resistance heating strips were 
cemented In a spanwlse direction spaced 1/32 Inch apart. A cowering 
of the l/61»— Inch-thick plastic— ln5»regnated sheet was laid over the 
resistance strips, and on top of this was cemented a skin of 
0.006-lnch-thlck aluminum. Each 1/2— Inch-wide heating strip was 
connected to individual liigs located along the edges of the test 
section. This provided means for chordwise adjustjaant of the power 
disliribution by l/2-inch increments. The heated area of the test 
section extended back to 77 percent chord on the left side and to 
17 n ercen t chord on the right slde^. 


The guard sections were constructed in the same manner as the 
test section, with the exception that the ftluTnimirn alfln was 0.011 
inch thick. The heated area of the guard sectione extended to I7 
percent chord on both sides of the model. 

Msasurements of the temperature of the umimini surface of the 
test Section were obtained by means of thermocouples. Fine iroi>- 
constantan thermocouple wire was rolled flat to produce a strip 
approximately 0.002 inch thick and I/16 inch wide. These strip 
thermocouples were laid in spanwlse grooves about 3 inches long cut in 
the aluminum skin. The thermocouple Jimctlons were located in the 
middle of the grooves, and the leads passed throu^ holes at the ends 
of the grooves Into the interior of the model. Aliamlnum was sprayed 
into the groove over the strip thermocouple for a distance of about 
3/16 inch on either side of the Junction. Thus, the thermocouple 
Junction was bonded to the aluminum skin, allowing accurate surface- 
temperature aseasurements to be made. The remainder of the groove on 
either side of the aluminum spray was filled with a nonelectrically 
conducting material, miermocouples were located at the center of 
the test section at 1— inch chordwise intervals In the leading-edge 
and calculated transition regions, and at 1-1/2— inch chordwise 
intervals In other regions. Surface ten5)eratures were recorded by 
means of self-balancing automatic— recording potentlometere. 

The flow of heat throu£^ the outer surface was calculated from 
measurements of the power dissipated In the electrical heating 
strips. This power was determined by meas\n*l33g the resistance of 
the strips and the current flowing through them. Thermocouples 
placed on both surfaces of the plastic base at a number of chordwise 
stations gave an indication of the heat flow into the mo del interior. 
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These thermocouples vere connected to the same recording potenti- 
ometers tised to record stnrface tenqperatvires. 

Pressure taps were installed flush in the test-sectioi surface 
about 3 inches down from the top edge of the test section at various 
chordvise points for the purpose of measuring surface pressure dis- 
tribution* A standard HACA 60 -<;ell pzassure recorder was used to 
record the pressures. 

A source of 400 -cycle, single— phase ^ alternating current was 
supplied to the test and guard sections for heating these surfaces* 
The heating strips for the test section were grouped into 30 chord- 
wise blocks* Control of the current flowing through each block was 
provided, so that a large variation in the chordvise distribution of 
heat flow was possible during flight* Before the icing operations 
stairted the heating strips for the guard sections were connected to 
maintain a constant sxa^ace— temperature rise during flight in clear 
air* Controls were provided so that the total heat input to the guard 
sections, but not the chordwlse distribution, could be varied during 
flight* 


HACA 0012 Airfoil Model 

With a few exceptions, the constiructlon of the 0012 airfoil 
model was substantially the same as that of the 6^,2-016 model* These 
exceptions will be noted* 

The top layer of plastic— la5)regnated fabric covering the electri- 
cal resistance stripe constituted the outer skin of the test and guard 
sections* This was painted and sanded, after the test-section thermo- 
couples had been Installed. The heated area of the test section 
extended back to 58 percent chord on the left side and to 11 percent 
chord on the ri^t side * The heated area of the guard sections 
extended to 11 percent chord on both sides of the airfoil. 

Strip thermocouples of the same type as Installed in the 65,2-016 
model were used to measure surface teu^^ratures of the tost section* 
Spanwlse gprooves were cut in the plastic— In^egpiated fabric sheet at 
various Intervals along the chord* The strip thermocouples were laid 
in the gprooves and cemented in place* The sinrface was then xn.inted 
and sanded so that only a thin layer of paint covered the thermo- 
couple Junctions* Thermocouples were located at the center of the test 
section at 1 -lnch chordwlse Intervals in the leadlng^dge and cedcula- 
ted transition regions, and at 2- to 2 - 1 / 2 -lnch chordwlse Intervale in 
other regions* Surface teaqaeratin^es wore reccn^ed by means of an 
automatic-recording potentiometer* 

<■ 

The flow of heat through the surface of the test section was cal- 
culated from measurements of the power dissipated in the electrical 
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heating strips. In a manner slmlleu* to that described for the 6^,2-016 
airfoil. 

Installation of surface pressrxre taps for the msasis^ment of ^ 
pressure distribution was the same as for the 6^,2-016 model. Ibe 
pressures were recorded by photographing a multiple— tube manometer 
to which the presstire taps were connected. 

A source of UOO-cycle, slngle-|hase, alternating current was 
supplied to the test azid guard sections for heating these sinr^aces. 

In the test section, provisions were made for obtaining a limited 
number of chordwise heating distributions as well as for control of 
the total heat Input with each distribution. A small degree of varia- 
tion of each heat distribution was also provided. During flight It 
was possible to control only the total heat Input, and to vary, to a 
small extent, each distribution. As with the 6^,2—016 airfoil, the 
-heating-strlps-for- the- guard^sejctions^jTO^e conne<^d to ^iye an 
approximately constant surface-teniperatin*e rise in cleair air. No 
control of the heat distribution or the total power Input to the guard 
sections was provided. 


TEST PROCEDIEE 

The test airplane was flown Into natural— Icing conditions over 
most of the northwestern area of the Ihiited States during the winter 
of 1945—46. During the winter of 1946-47 the area of operations was 
extended to Include a few flints In the central and eastern part of 
the Ihilted States. The usual test procedure, during fli^t in Icing 
conditions, was to records airfoil data simultaneously with the 
measurement of the meteorological conditions. The rotating cylinders, 
described In reference I3, which ccxistituted the means of measuring 
llquid-^ter concentration, drop size, and drop-size distribution, 
were extended as often as was conveniently possible. Becords of 
free-air temperature, airspeed, and altitude wars taken sevez*al times 
a. minute. The recording potentiometer xised to obtain airfoil tenpersf- 
tures was operated continuously. During this time, the values of 
current flow throu^^ the electrical heating strips of the alz^oll 
were recorded. Fhotographe of the test-section surface and records 
of pressing distribution were taken at freq,xient intervals. 


RESULTS 

A tabulation of the fli^t and meteorological conditions for 
which simultaneous airfoil data were obtained Is presented In tables I 
and II. Table I contains the fli^t and icing conditions for which 
corresponding heat— transfer measurements were made with the MCA 0012 
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airfoil. Table II gives similar information for the NACA 65,2-016 
airfoil. All measurements vere made during fli^t in nature^.— icing 
conditions. During most of the flights, large variations in liquid— 
vater concentration, and occasionally drop size, were experienced. 

The rotating cylinders, used to measure liquid-water content and drop 
size, were extended for about 1 to 2 minutes, thus giving average 
values for 1— to 2-mlnute intervals. A conq)lete cycle of the 
recording potentiometere used to record airfoil tenqjeratures required 
h to 6 minutes, during which time the meteorological conditions may 
have changed considerably. For these reasons, an effort was made to 
select, for analysis and discussion herein, only the airfoil data 
recorded during flight in relatively uniform clouds and/or \diere 
close correlation existed between the cylinder measurements the 
airfoil— tenq)erature records. Of these data, only a part, chosen as 
being typical, are presented in this report. These are the thermal 
data for fdilch the fll£^t and Iclzig conditions are given in tables I 
and II. 


NACA 0012 Airfoil Data 

Figures 8(a) to 8(g), inclusive, present the measurements of 
siirface tenq>erature, surface heat flow, and resulting heat-transfer 
coefficients obtained with the 0012 airfoil model during fll^t in 
the conditions presented in table I. The heat-flow distribution 
Illustrated in these figures had been found by experiment to give 
an approximately uniform teBq)erature rise over the test-section 
surface during fll^t in clear air. Variations in the Intensity of 
the distribution for the different conditions of table I occurred 
as a resTilt of the heat supply procedure followed during the teste. 
In general, during an icing test the total heat input was reduced 
tuitil the surface temperature was observed to fall close to freezing 
tenqperature at some point on the test section. Typical values of 
eiurface temperature, heat flow, and convective heat-transfer coef- 
ficient obtained during fli^t in clear air are shown in figure 8(h). 

Qhe data presented in all figures except figure 8(g) were taken 
with the entire test section heated. Figure 8(g) presents data 
secxired with only the leading-edge region heated, from 11 percent 
chord on the ri^t side to 8 percent chord on the left side. At the 
time of this test, insufficient heat was supplied to the leading- 
edge area to evaporate all the water striking the surface, and 
streamers of ice formed £ift of the heated region, similar to those 
shown in figure 9. 

The heat— flow values given in figures 8(a) to 8(h), inclusive, 
were calcvilated from measurements of the totsil power dissipated in 
the electrical resistance strips and the internal heat loss. The 
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meastEPeinents of surface tell^)erature4 idilch vere oT)tained with the 
thermocouple Installatlotu prevlo\isly described, were corrected for 
errors Incurred by the presence of a layer of paint covering the 
function. Magnitude of Jbhese enrors was determined from knowledge 
of the thickness and thermal conductivity of the paint and the amount 
of heat flowing throu^ lt« 

The kinetic ten5>erature of the free-stream air t©^ used In 
conqouting the heat^tranefer coefficients, was calculated for the 0012 
airfoil from equations (I8) and (I9) , using experlmenteOly determined 
values of the espresslon 


l-^(l-ErP) 

— ^where — 3 — Is- 1 / 2 - f or-lamlnar-flow-and-l/3 f or— turbulent-flow.-Value s- 
of this eapresslcn for various points along the airfoil surface, were 
obtained from figure 10 , which presents data obtained during fll^t 
In clear air. 

A typical record of pressure distribution over the 0012 airfoil 
model test-section surfeuse Is shown In figure 11 , 


HACA 65,2-016 Airfoil Data 

Figures 12 (a) to 12 (j), Inclusive, present the measurements of 
surface tenq)erature, surface heat flow, and resulting coefficients 
of heat transfer obtained with the 65,2-016 airfoil model during 
fll^t In the conditions presented In table II. The distribution 
of heat flow shown In these figures had been e^>erimentally estab- 
lished to provide an approximately uniform tenqperature rise above 
free-alr ten5>erature over the test-section surface during fll^t in 
clear air at an altitude of 11,000 feet and a true airspeed of I75 
miles per hour. This heat-flow profile was used throu^out all the 
fll£^t tests. Slight variations in heating intensity are due to 
variations In Internal heat flow and chordwise heat conduction in 
the thin altaninum skin. Typical values of surface temperature, 
heat flow, and convective heat-transfer coefficient for flle^t In 
clear air are given In figure 12 (k). 

The results shown In all figures (except figs. 12 (h) to 12 (j), 
Inclusive) were obtained with the tost section heated to approximately 
55 percent chord. Heating aft of this point was precluded by a 
malfunctioning of the heating equipment In this area. Figures 12 (h) 
to 12 (j), Inclusive, present data secured with only the leading-edge 
reglcm heated, frcmi I7 percent chord on the ri^t side to I7 percent 
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chord on the left side. During the time vhen the measurements of 
figures 12(1 ) 12 (j) vere taken^ an insvifflclent quantity of heat 

vas being svqpplled to the leading-edge region to evaporate all the 
vater striking the surface, and ice accumulated aft of the heated 
area. This is shown in figure 13* 

The heat— flow values given in figures 12(a) to 12 (k), inclusive, 
vere calculated from measurements of “^e total power dissipated in 
the electrical resistance strips, in a similar manner to that used 
for establishing the heat flow for the NACA 0012 airfoil. In 
addition to the determination of the internal heat loss in coioputlng 
the surface heat flow for the NACA 65,2— 0l6 airfoil, the flow of 
heat chordvlse in the thin aluminum surface was considered. The 
chordwise heat conduction is a function of the chordwise variation 
in surface temperature. It was assumed that the quantity of heat 
flowing from point to point eilong the surface, as indicated by the 
difference in surface temgperature between the two points, originated 
from the heating strip under the hig^r tenq)erature and flowed away 
from the sutrface into the air stream in the area of surface over the 
heating strip at the lower temperature. This method, althou^ 
inexact, offered a rapid means of estimating the effect of chordwise 
conductlcui. A more exact determination of this effect can be 
obtained using the ’'relaxation" method of reference 10. Ho correc- 
tions were applied to the su3rface— teoqperature measurements, since 
it was assumed that the sxxrface thermocouple Junctions were at 
surface tenq)erature. 

The kinetic teuq^rature of the free-stream air to^ was 
calculated using equations (I8) and. (19)* Values of the e:q>ression 
1 - (l?/V®) (1-ErP) were calculated and are plotted in figure l4. . 

A picture of the conditions of wetness ^ich existed on the 
airfoil during flight in clouds can be seen in figure 15. This 
figure shows some typical records obtained with strips of blueprint 
paper which had been fastened to a device that could be extended 
into the air stream up the leading edge of the airfoil model to a 
point Just below the test section. Since, in effect, these vere 
wrapped around the leading edge of the model, they illustrate the 
pattern that the water ass\anes in striking the airfoil and flowing 
aft. The records were obtained during icing conditions 11, I3, and 
14, table II. 

A typical measurement of pressure distribution over the 65,2—016 
airfoil model test— section surface is shown in figure I6. 
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DISCUSSION 

I3ie ultimate in performance of a vlng thermal ice— preTentlon 
system Is one idilch vUl prevent the accretion of ice cm. any portion 
of the vlng* Qhls Ideal operation requires that any vater on the 
vlng surface must he maintained in a liquid state until it evaporates, 
or hlovs off the vlng at the trailing edge* In many vlng designs, 
heating of the entire surface is not practicable because of such 
features as integral fuel tanks, and in these instances any vater 
flovlng aft of the heated region is apt to Areeze wand fcarm the type 
of ice accretion normally termed runback* In the follOvlng dis- 
cussion the ice-prevention action of a heated vlng vlU be examined 
in detail, and the reliability of the equations assumptions 
presented in the analysis sectl^ for the prediction of surface 
temperature and rate of vater evaporation from the surw.ace vill be 
established* If these equations and emalytlcal methods can be shown 
to 'deflne~correctly-tho-procoss-of--thermal^lce_preventiOT, "toe fundsi- 
mental design procedure for a heated wing initially ctaiceived in / 
reference 6 will be more firmly established* Ihe en^lrlcal design . 
method of providing a specified ten^rature rise in clear air can 
then be replaced by the more fundamental and flexible concept of 
supplying sufficient heat to maintain the surface temperature above 
freezing until the vater is either evaporated or carried away* 

An emaLlysls of the action of a heated vlng requires the consid- 
eration of three factors: namely (l) the meteorological and flight 
conditions for idilch the vlng must provide protection; (2) the area 
of vater in^ingementi and the rate and distribution of li^plngement 
over that area; and (3) the rate at vdilch the vater is evaporated 
from the vlng surface* 


Meteorological and Fli^t Conditions 

The specification of a meteorological condition for the design 
of thermal ice— prevention equipment depends upon the geographical 
areas over vdilch the airplane vill fly, the seas(ms of operation, and 
other factors dictated by the intended service of the aircraft* 
Obviously, the establishment of design conditioxis for a specific area 
requires a knowledge of the conditions prevailing over the area* If, 
on the other hand, the ice— prevention system is to provide protection 
for all-weather operation, general specifications of a meteorological 
condition must be established which vill encoBq>ass all conditions 
likely to be encountered* 

Ihe most recent and extensive informaticui in regard to the 
severity of icing conditions likely to be e:^rleiu:ed in alL-weather 
deration in the IMited States is contedned in references I3 end 14* 



22 


HACA TN Ko. 1472 


In reference 13^ estlmteB of the wt^Tiwuin contlnuone icing condi- 
tions as veil eu 9 the maxlimna prohahle icing conditions apt to he 
encountered are presented. Since the duration of the 
prohahle icing condition is quite short (1 to 2 oinutes)^ and, icing 
of this seyeritj is entirely associated vith cunuliis clouds vhich 
should he^ avoided in all operations^ the aftTimmi .p-nrit iniifH ^a icing 
condition' is helieved to ho of greater interest for design purposes. 
Tvo conditians of maxlimiin continuous icing are presented based on a 
relationship of drop size and liquid-^iater content. Ihese conditions 
are given in the following table: 


Liquid-water 

concentration 

(©a/m®) 

^Mean-effective 
drop diameter 
(microns) 

Free-air 

tenperaturo 

(°p) 

0.8 


20 

0.5 

25 

20 


It is helieved that the conditions in the above table form a good 
basis for the design of thermal ice— prevention eqtdpment for «n— 
weather operation. In addition to these values, however, the proposed 
wing thernBl system should he analyzed for possible undesirable oper- 
ation in other icing conditions. For example, reference 13 points out 
that drops of 35 to 50 microns diameter should not be regarded as 
exceptional. Althou^ the amotmt of liquid water associated with such 
large drops is usually low (about 0.1 gram per cubic meter) the fact 
z*emalns that the area of water liiq>ingenent would be very large 
would probably exceed the limits of the heated region if this region 
had been based only on a consideration of the data in the waTiTmitn 
continuous table. F ina lly, the possibility of encoimterlng icing con- 
ditions at low tenperatures nay be a critical condition for heated 
wings on same airplanes. For Instance, the estimated conditions 
of naxlmnn continuous icing presented in reference I 3 ana given 
in the preceding table were extended in reference l4 to air teuqwra- 
buTOs as low as -50® F for the case of 15—micron drops. The coit- 
dltlons of maximum continuous icing suggested are 0.5 gram per cubic 
meter at 0° F and 0.25 gram per cubic meter at -50® F, both with 

Ihe msaiv-effective diameter as defined in reference I 3 is the size of 
drop in a cloud sanq)lo for which the amount of liquid water existing 
in water drops larger than that drop is equal to the amount of liquid 
water existing in drops smaller than the drop. 
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a maan-effecti-ro drop dlaneter of 15 microns. Either of these condi- 
tions may he more deleterious to the functioning of the ving thenoal 
system than those in the table. 

From the foregoing discussion^ it is evident that the analysis 
of a heated wing should give consideration to Several possibly criti— 
ce^. icing conditions in the same manner that severed, flight condi- 
tions are assumed for the ving structured, analysis. The data of 
references I3 14^ althou^ somevhat limited in scope, are consid- 
ered to be sufficiently Indicative of icing coxiditlons in the Ihilted 
States to form a meteorologlceil basis for heated-^i]g design. 

Thft problem of selecting a f li^t condition for the design of 
ice-prevention equipment is concerned vith the airspeed and altitude 
at ^diich the airplane vlU fly. The airspeed will depend upon the 
specific airplane, and, in general, a cruise condition should be 
“selected^ — Chooslng an altitude-fQr-deslgn-ie-dependent_upon_ee.veral_ 
factors, vhlch will be discussed later. 


Area, Bate, and Distribution of Water Iiaplngement 

Saving defined the icing conditions for idiich the heated ving is 
to be designed, the next step is to determine the region of the 
leading edge in idilch the water drops will strike the wing, the rate 
of water iziq>act at any specified point in that region, and the total 
rate of impingement per foot of ving span. This subject was dis- 
cussed at some length for the general case in the analysis section. 

In that discussion, it was shown that the method of reference 15 could 
be used to prepare (for any wing section for ^dilch the stream lines 
wez*e known or could be determined) curves similar to those presented 
in figures 1, 2, and 3. The broken line of figure 1 gives an indlcai- 
tlon of the area of water isplngement, vdille the rate of Impingement 
at a specified point can be obtained from figure 2 and equation (7)« 

Two methods are available for the determination of the total 
rate of Impingement per foot of wing span. The calculation of , this 
quantity is of primary importance, as it determines the amount of 
heat required to disperse the water by evaporation. The first method 
utilizes the concept of collection efficiency E as mentioned in 
the analysis section. This method is preferable vhen only the value 
of the total rate of impingement is desired, since preparation of 
the curves of figure 2 is not required. For a thorou£^ analysis of 
the heat transfer from the surface, however, knowledge of the rate 
of iaplngemant at a point % is required. By enploying equation (7) 
and figure 2 a curve of the distribution of water inplngement 
(1% against s/c) caii be plotted. Figure I7 shows such a curve for 
the NACA 0012 airfoil, using eqmtlon (8) and an ”E" type drop-size 
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distribution. (See reference 17.) A curve of this type presents an 
interestl33g picture of the distribution of vater liiqplngeiaant and^ in 
addition^ the area under the curve denotes the total rate of vater 
interception. 

Althouf^ the method of reference is considered to provide a 
coD9>lete and q.ulte accurate prediction of the distribution of vater 
la^lngement on the leading edge of an airfoil^ it does have the dis- 
advantages of requiring (1) a knowledge of the velocity con^jonents 
along a nuoiber of the airfoil stream line s^ and (2) considerable 
computation. The difficulties associated with the cooqiutation of 
the vatei>-drop trajectories for airfoils have encouraged the substl— . 
tut ion of a cylinder vlth radius equal to the airfoil leading-edge 
radius in the determination of vater impingement. (See references ^ 
and 6.) Ohe curves of reference I7, \dilch have been calculated for 
a large range of drop sizes, airspeeds, altitudes, and cylinder 
diameters, are then used directly to evaluate the anticipated vater 
inqilngement on the airfoil. This substitution procedure is a useful 
device but should be esployed with a full knowledge of its llmitar- 
tions. One of these limitations is the fact that the curves of 
areference I7 provide the area and total rate of water iB^plngement, 
but give no direct indication of the distribution of lmplngea»nt. 

A second restriction of the cylinder-substitution method is 
concerned with the contour and size of the forward portion of the 
airfoil. To obtain an indication of this effect, the rate and area of 
vater lnqplngement on the 0012 airfoll,^ at 0° angle of attack, and on 
the leading-edge cylinder of that airfoil are conq>ared for the same 
fll^^t ccmditlons and various drop sizes in figure I8, Ihe values for 
the 0012 section were obtained from flguzss 1 AT>d 2, and those for the 
cylinder from reference I7. At drop diameters up to about 25 microns the 
rates of is^ingement on the airfoil and on the leading-edge cylinder 
are approximately the same, althotj^ above 25 microns as the drop size 
Increases, the rate of isq>ingement becomes considerably greater on 
the aiz*foll than on the cylinder.. At drop diameters up to about I8 
microns the area of in5>lngement on the cylinder is roughly equal to 
that on the airfoil. However, at a drop diameter of 25 microns, which 
is not musual (reference I3), and was presented previously in this 
report as a possible maximum continuous condition, the area of inplnge- 
ment on the airfoil is neau*ly 50 percent greater than on the cylinder. 
It should be noted that the value of 25 microns for the nw-Y^nmni contin- 
uous condition is the mean-effective diameter, and that drops of a 
larger size probably will be present duo to the existence of a distri- 
bution of sizes. Although these valtuas provide an Indication of the 
scale limitation of the cyllnder^ubstltution method, the fact should 
bo noted that figure I8 applies to only one airfoil section, with 
an 8-foot chord, and at one fli^t condition. The leading-edge radius 
of the NACA 0012 section for an 8-foot chord is small (1.5 in.) 
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and the leadlng-ed^ cylinder does not oatch the section contour for 
any great extent above the chord line. This Is shown graphically In 
figure 19 which presents a coBqnrlson of the forward portions of 
three airfoil sections and the leading-edge cylinder of the 0012 
section. In the case of edrfoU sections with the leading-edge radius 
a greater percent of the chord than the 0012 section, also for 
airfoils of 0012 section, or similar, with chords greater than 8 feet, 
the cylinder>-substltutlon method will present a better approximation 
than that Indicated by figure I8 for the sane speed range. 

For airfoil sections with a leading-edge radltis idilch represents 
a small percentage of the chord, the substitution of an "equivalent" 
cylinder (referezuse 12) with a radius larger than the leading-edge 
radius would probably provide a better Indication of the rate and 
area of water laqplngement on the airfoil than would be obtained for 
the leading-edge cylinder. At the present tine there Is not suffic- 
ient Izjformat Ion (mwater-drOp-tra^ectories about-alrfolls_to_provlde_ 
a basis for selecting the proper cylinder In each instance; therefore, 
the designer must utilize the more coaqilicated, but more accurate, 
method of reference 1^ <3r assume some cylinder diameter based on his 
experience. The possibility that the rate and area of laqalngenent on 
an ellipse would more closely approximate the rate area of 
Impingement on a series of similar airfoils has been suggested 
Is worthy of future consideration. 

The ability to select a proper drop size for the design of wing 
Ice-^preventlon equipment Is a factor of considerable Importance to 
the designer, as can be Illustrated by figure 20. In this figure the 
rate and area of lnq^lngement are presented for the 0012 airfoil as 
a function of drop Size. The rate of laplngement for each drop size 
was ceiLculated for a llquld-^iater content of 1.0 gram per cubic meter. 
Consider, then, a change In design drop diameter fjrom 10 microns to 
20 microns. The resultant Increase In rate of water iiqilngement is 
1.7^ poxmds per hour per foot oi span or an Increase of 17^ percent, 
al-^ough the actual amount of water present per xmlt yolume of cloud 
has not been changed at all. The same Increase In drop size will 
cause an Increase in area of impingement frcm 1.^ to k percent s/c. 

In contrast, consider the effect on the rate of water Impinge- 
ment produced by an Increase In the quantity of liquid water present, 
assuming the drop size to remedn constant. The area of Isqtlngement 
will remain unchanged^ idiUe the rate of lsq>lng«nent will Increase 
only In direct proportion to the Increase In water concentration. 

This ex82i^)le clearly Illustrates, the fact that the amount of free 
water p^sent In an icing cloud is only ox» factor Influencing the 
quantity of water idilch will actually strlhe the wing In a specified 
time Interval, that the size of the cloxad drops Is a factor of 
at least equal inq>ortance. 
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The prohldu bf dlstrihutibh of tkb slzbs of dfppe ^ ah. icirig 
cloud alsb hecirs careful cb^idefaiioh^ ^or eiliii^iei if ah typb 
drop-elze dlatrihutibh exists a hebni=effbcii^ drop dianhter 6f 
2^ Bdcrons^ largest drbps vlli hb ^ ^crbhs in diameter* ibhce^ 
the area of itipizi|^mBht id.ll he cbhsiderahly heater t]^ if a 
uniform drop size of 0 microns prei^ied* The data bf reference l4 
indicate that> in general, the dlstrihutibhs of drbp size in icing 
clouds are fairly haarbv> and do not udually f ollov the broad distrl- 
butlons, such hs type Nevertheless, the dlstributicai Dhiet be 

considered^ since the lar^st drOps in the cloud determine ti^ area 
Of impingement and the mihimim extent Of heated area required for 
ice preTentlOni 


Rate Of ErepOwition of Water 

Having discussed the problem of ai^a and rate of Water intercep-; 
tion, the next step is to establish the rate at idiich the intercepted 
water is evaporated from the vihg surface and the vaiidity of the 
eq^uations presented previously for determining t^ rate of heat 
dissipation during the process of evaporatichi Tbe problem of rate 
of evaporation is partibuiariy impcirtant because all Of the Water 
intercepted by a wing heated Only- in the re^oh Of the leading ei^e 
must be dispersed by evaporation if the fOrmtion of rUnbach is to 
be avoided; 

From a superflciai study Of the mechanism by Which Water is 
deposited on the surface of a wing; it WOiild be eapected that in the 
area Of water iapingement the surface is cOnpletely Wettedi and that 
equation (2h) for calcuiating the hCat loss from a heated wing is 
valid. Aft of the area of in^ingement; it Would be anticipated 
that the surface inay not be fully wettedj Since water dOes not reach 
this region directly^ but instead must flow bach ffcm the area of 
impin^ment* If the Surface aft Of the area bf impin^ment is Only 
partially wetted, the eacpression fin* i (eqUatibn (0)) iffiist be 
modified for use in equation (25) tO calculate heat requirements • 

Observatlons> made during the airfoil tests. Of the water 
pattern On the leading ed^ Of thO alrfOils revealed that the above 
ex^posltloUs are Correct; At a very short distance bach of the 
re^^On of Inpin^ment, the film of Water was observed to reach a 
state Of lilStabillty and break iUtO small rivulets; A picture Of 
the con^tibnS Of Wetness which actually exist bn a Wiiig during 
flight in clOUds Can be seen in f ig^ 15» , ^is figure, which sHws 
typic^ recOTds Obtained With the strips Of blueprint paper placed 
around the lee4ing ed^ Of the 0i2-iJDl6 airfoil model during fll^t 
In icing conditions, illustrates the patterii fOi^d by the water in 
striking the airfoil leading edge and flbwing aft. It i§ evident 
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from the patterns that the area of ln5)lngeinent, which Is clearly 
defined. Is con5>letely wetted, while hack of this area the water 
collects and forms rivulets, creating a partially wetted surf ewe* A 
study of the patterns indicated a variation in the fraction of surface 
area wetted (aft of the impingement area) with rate of impingement of 
water. Accordingly, the rates of water impingement Mg were cedcu— 
lated using equation (U) for the conditions existing at the time 
that rivulet patterns were obtained (icing conditions 11 throvi^ 15, 
table H). The curves and values presented previously for the OOiP 
airfoil were used in the calculations of Mg for the 65,2— 0 l 6 airfoil. 
Substitution of the calculations for the 0012 section in conputing 
values for the 65,2— 016 section appears to be a good approximiatlon, 
since, the contour of the 65,2—016 section in the leading-edge region 
is very nearly the same as that of the Joukowski airfoil used in the 
0012 trajectory calculations. Figure I9 compares the contours of 
_the_three sections. 

The values of Mg were plotted against the measured areas of 
surface wetted, obtained from the strips of blueprint paper. 

Figure 21 shows the relationship, thus obtained, between the rate of 
flow of water from the impingement region and the fraction of sttrface 
ai^ea wetted. For the data shown in figure 21 , values of the rate of 
water flow over the siarface aft of impingement were assumied 
to be equal to the rate of water Impingement Mg, Ihe scatter of data 
points in the figure is believed to be caused by errors in measure- 
ment of the liquid- water concentration occurring at the time the 
rivulet patterns were obtained. Table H shows that the free-elr 
temperature was hi^ during icing conditions 11 throuf^ 15, •»dien the 
blxieprint records were taken. The kinetic temperature was close to 
freezing, and it was observed that the water striking the rotating 
cylinders, used in the measinrement of water concentration, was 
running back, and possibly off, the cylinders. Thus, the liquid- 
water concentrations measured may have been lower than the actual 
concentratione present, Ihe two data points corresponding to a 
wel^t rate of water flow of 0,57 pounds per hour per foot of span 
(fig. 21) irepresent the rivulet data procured at the lowest free- 
air temperat\ire of these tests (icing condition 14 , table II), Ihese 
are probably the most reliable data, since the rotating cylinders 
were subject to smaller losses of water, Kierefore, the curve shown 
in figure 21 was wel^ted toward these points. The ultimate extent 
of this curve in the direction of percent of surface wetted is not 
definitely known. Ihere is evidence, however, indicating that the 
degree of surface wetness aft of the area of impingement reaches a 
maxlmnim which is not exceeded, regardless of the rate at which wa'tor 
is Intercepted. It was found that a relationship exists between the 
rate of flow of water in the region aft of impingement and the 
surface-temperature rise above free-air temperatiu?e, and that the 
temperature rise decreases to a llmilt as the rate of flow of water 
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Incre^es. Figj^e 22 shows the ^lati<^hip l^tmen the rate at 
which flo^ hi^k over Idle heated syrf^e of the 65,2:^16 aiiw 

foil test eectlon and the syerage i^ficase in teiqperature ahove 
fre.e^if pf the surface ^oin 10 to 2^ ^rcent chord. Ihe 

T^j^s pf the ws'ter ^ow were obtained hy suhtractipg from 

tfci calculated rates of rater implngenent for one side of the airfoil 
the corignatod raH»s of evaporation from the region of ln^lngement^ 

Figinre ^ illustrates that as the rate of watjer flawing over the 
sut^ace incraaseSj t^ of 1^ surface decraases, hut that 

a llMt to t^ decrease in tei^ratiure appiu*en^^ is reached. iSils 
Indica'tos ti^t rate pf e'vapQration reaches a sEucimum as the 
limitl^ surface tei>S!eratlih^e is apprpechedj, since evap<n‘ation is 
the pi^ veriahie in the h©at--transfer process in the area of surface 
uMer "cpnsitoration. -Therefc^, if the ra^ of eyappratipn attains 
a muiimm degree Pf surface wetness must ^sp approach a limit, 
it cw 1 m dem^^ using equation (25) ai^ i^e v^ues presented 

in f Igi^s ^ aiid 22^ that the maximum fraction of surface area 
ratted is ahput 5P ppraent. 

Mthou^ the 4^^ frcm which the curve of figure ?1 ras confuted 
rare phtalnSd with blueprint paper strips wrapped arpu^ the leading 
edge pi 65^2-^16 airfp^ giyen in this figure are 

■believed tp 'be sufficiently indicative pf'^ the c, editions of wettability 
existing ph all clean wing surfaces not specially traated to be appli- 
cable for general a^^Pti tl^rmal desij^«r^ For purppses pf design^ 
it is sugessted that the lii^t 9f surface ratneSs foy surfaces not 
special^ ■fciea'ted 'be taken as hp ^raeiit* It is pf in^crtance to 
npte thai' "the Cb?ve pf surf ape ratness shown in figure 21, 

for a heated :^ng, the total ra'te of eyappratipn of rater Ve in 
the reglm pf raterr^op Ir^iiiigeiaent must subtracted from' the 
totai rate Pf rate? implngemen't in order to obtain the rate of 
flow of rater rearrard from the area pf impingement. ' The values 
giyen in fign’e 2I for degree pf surface l^tness'are 'beileyed to 'be 
accurate phly'tp tho lissrest IQ percent. 

With the informatipn gained sp far, it shpuld "bp ppeslble to 
analyze the da'!^ obtained with the tVP electrically heated airfpil 
models ai^ esta'blish tixp y^ldlty pf equations' (24) and (25) for 
calculating ^at flow. The curves of masured heat-flow distribution 
shown in figures 8(a) tp 8(g) flgi^s 12(a) to i§(j) rare faired 
to prpduce e form more sul'^'ble fpr bpnq^ispnwith hea"!^ curves 
c^cuia^d tiBing equations (24) and (25). Qp^arispns of the 
measured heat flow and the heat flow peculated to ^oduce the 
measxirad surface i»nperatures, assuming tirn'entlre surface. tp be 
cpj^leitely rati^d# fbr the t'^o airfpil sections fey t^ical cases eye 
Shcra in figures and ^4. These cu^S srp^^’s^ cpjf^^eiii tp the 
CJ^CUlafed heat los due tP PPUTPoti^ pi^y, is, assuming the 
s^face to be cp^ietely dry (equatlOT (ll'))» Ip ; the previously 
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msntlofned calcttlatloos, wasured Taluee of convoetlve healv-transfer 
coefficient, obtained during flight in clear air, vere used. In 
order to calctilate the amount of heat dissipated in warming the 
Impinging water (equation (4)) for the calculation of heat flow for 
the completely wetted case, values of were conqputed from equar* 
tion (8) ttsing the measured values liquid-water concentration, 
drop size, and droi>-elze distribution (tables I and U). As was 
done previously in the analysis of the blueprint-paper rivulet 
patterns, values of % were computed for both airfoil sections 
using the curves presented for the 0012 airfoil. ' 

A study of the measured and calctilated heat-flow curves In 
figures 23 and 24 shows that In the area of water-drop li^lngmaent 
good agreement Is obtained between measured values and the values 
calculated for a completely wetted siarface. Indicating that In the 
reg ion idiere It Is reasonable to fissume a fully wetted surface the 
equatlOTs for calculatihg^heat^flow are vall — Aft-of-the-area_ o f 
Impingement, In the region of low heat flow, where It has been shown 
that the surface Is only pasrtlally wetted, the values calculated for 
a completely wetted surface are lower than the smasTxred values. Since 
the BvarSace Is only partledly wetted. It would be eipected that the 
calculated curve, idilch represents the values of heat flow required 
to produce the measured surface temperatures If the su^ace were 
completely wetted, would be considerably higher 'ttian the measured 
curve. There appear to be only two possible eaplanatlons for this 
discrepancy: (1) equation (2^) gives erroneous values cannot be 
relied tpon for calculation, and (2) the values of convective heat- 
transfer coefficient used In equation (2^) for calculatix^g the heat- 
flow venues are In error. Ihe first explanation does not appear to 
be likely In view of the fact that the equation was derived on a 
sound basis. Also, there is no obvlonS reason lAy the equation 
should hold In the leading-edge region and fall to hold in the area 
aft of the leading edge. Ihe second esQtlanatlon, that erroneous 
values of convective heat-transfer coefficient were used In the cal^ 
culatlons, seems entirely possible. Since the valtms of convective 
heat-transfer coefficient used In equation (2^) are those measured 
during fU^t In clear air. It seems reasonable to assume that tran- 
sltlcm ft*om laminar to turbulent flow moved forward during flight In 
Icing conditions frcai the position maintained In clear air. Movement 
of transition to a point near the leading edge would cause the ccuxveo— 
tlve heat-transfer coefficients In the region under consideration to . 
be Increased several times above the values existing In clear air, 
since the convective coefficients In turbulent flow generally are 
considerably greater than those In laminar flow. Such cui Increase 
In the convective coefficients would raise the curves calculated for 
a completely wetted surface (flgs« 23 and 24) to a position above the 
measured, curves. / 
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It is ■ very lllcely that dlsturhazice of the ho\mdai7' layer, caused 
hy water drops striking the airfoil surface and roughening the surface 
as they coalesce and flow aft, would effect a forward movemant of 
transition. Ihere is further evidence to support the assunqptlon that 
the watezv-roi^ened surface catised movement of transition forward. 
Observations of the 65,2-016 airfoil during a fll^t in clear air 
with the test section heated Indicated that transition hflVi shifted 
forward by a considerable amovmt. This was noted by a lowering of 
the surface tenqperatures in the region aft of the leading edge. The 
' heat distribution had been set previously to produce a constant 
surface temperature in clear air, and only a change in the boundary— 
characteristics could cause the evident change in heat— transfer 
coefficient. After the fli^t, a close examination of the leading- 
edge region of the airfoil revealed small insects stuck to the 
surface where they had hit during the flight. The surface was wiped 
clean and during a subsequent flight in clear air it was noted that 
transition had moved hack again, as evidenced by the restoration of 
the surface temperatures to nonnal. ^us, it appears that very arnni ^ 
irregularities in the surface, such as are present on the sxnrface of 
an airfoil in icing conditions, are sufficient to cause transition to 
occur prematurely. Tests in wind tunnels also have shown that aTn«n 
protuberances in the leading-edge region of an aLirfoil will cause 
the movement of transition forward. (See reference 24.) 

Most of the curves of heat— transfer coefficient measured during 
flight in icing conditions, shown in figures 8(a) to 8(f) ani^ 12(a) 
to 12(g), display a definite increase in the aft region of hlg^ heat 
intensity, suggesting that transition is located at this point. It 
should be noted that the increase in heat— transfer coefficient 
indicated by these curves is believed to be only an apparent increase, 
caused by the rapid change in heating intensity in this region. If 
the coefficient is relatively constant throvi^out this area, as it is 
believed to be, a sudden increase in heating intensity will not be 
accompanied by an equally rapid cheinge in the thermal boimdary layer, 
and fcxr a short distance aft the indicated values of heat-transfer 
coefficient will be erroneously hl^. 

The exact values of the convective heat— transfer coefficient in 
the region aft of the area of iirq>ingenient in icing conditions etre 
unknown, but it is believed the values fluctuate due to changes in 
the location of transition during fli^t. Very probably, the 
disturbance to the boxaidary layer caused by water on the airfoil 
sxn*face is of such a character as to create instability in the 
boundary layer, and cause the location of transition to fluctuate. 

In the aft region of high heat flow (figs. 23 and 24), the 
values Of convective heat— transfer coefficient are Icnown, since 
tinrhulent flow existed in this region in clear air, when the values 
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vere neasirred, as well as in icing conditions. In this region, the 
measured heat— flow curve eind the calculated curve of convective heat 
transfer come together, indicating that at the point where the cvirves 
coincide all the water on the surface has been evaporated. 

Since the equations for calculating heat flow liave been shown to 
be valid in the area where the surface is canq)letely wetted, it is 
reasonable to assxxme that the equations hold in regions where the 
surface is only partly wetted, provided the correct modifications are 
made to the evaporative factor X. Fairly accurate modifications to 
the factor X are believed to be possible by using the curve of 
surface wetness shown in figure 21. By the use of this curve it 
should be possible to calculate the rate of evaporation of water from 
the surface aft of the region of water impingement. In the region of 
lii5)lngement the calculation of rate of evaporation is strai^itforward, 
since -full-evaporation_oo_cjurs»_ If the rates of evaporation from the 

two test airfoils can be demonstrated t^o be equiT^to' the Irates^of 

water inq)ingement for the test conditions, the method for calcula- 
ting rate of evaporation will be substantiated. 

Accordingly, calcxilations of the rates of evaporation from the 
surfaces of the two test airfoils were made for all the conditions 
of tables I and II for which thermal data were obtained. The rates 
of evaporation were determined graphically, using the curves of 
measured heat flow and calculated convective heat loss similar to 
those shown in figures 23 and 2k, Aft of the area of liiq)ingement, 
the position of the convective curve was established by dividing the 
measured values of heat flow by the modified values of X (equa- 
tion (25)). Values of the degree of surface wetness used in modify- 
ing X were determined from figure 21, using conq)uted rates of 
water impingement and evaporation from the area of iiiq)ingement. The 
position of the re— calculated curves of convective heat loss are 
shown typically in figures 23 and 2k, The total rate of evaporation, 
then, was determined by measuring the area between the mea.siu’ed 
re-calculated convective heat— flow curves. This gave the total 
amount of heat dissipated by evaporation of the water in Btu per 
hour per foot of span. Dividing this value by Lg the latent heat 
of vaporization, the total rate of evaporation Wg in poiuids per 
hour per foot of span was obtained. ' 

The rates- of evaporation, ob-tained in -the previously mentioned 
manner, are compared with the rates of water impingement, calculated 
by the method previously presented, for the 0012 and 65,2-016 air- 
foil models, for the left side only, in tables III eind IV. An 
average agreement of 13 percent for all the conditions analyzed 
where no runback formed was obtained, indicating the degree of 
reliability of the method for calculating the ra-be of e-vaporation 
of wa-ter from a heated wing. 
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In order to demoiwtrate fmrther the dependability of the ne^pd 
for calculfiiti^ rate of eyapprati^ from a heated wing, 
of rmbaph pn the 6%2~016 a^f^^^ I3) was ai^yzedV If it 

can be shown that the actual rate of formation of rimbacic c 
closely to the rate at lAich i^hack is caiciiiated to form under the 
particular icing conditions, tto meiihpd fear calci^ting rate of eyap- 
oration will be further subst^tlated. ~ . 

Ih® ?nuiback shown in figinre I3 had started forming 10 minutes 
earlier. At the time of the photograph the formation was estimted, 

"by obsei^tipn duripg fli^t, to he approrimately 3/I6 iiwh thipk. 

The area of the f oration exte;^d a^ut 2—1/2 inches cherdwise and 
12 inches spai^ise, making a wei^t'^of ice pf 0.2 poui4. This 
constitutes an ewtnai rate of foi^tion of runback of 1,2 pounds per 
hour per foot ep^. During this l6-minute peric^, tyP sets P? 
rptating-rcylinder and airfoil heat— transfer data were taken, ^eee 
correspond to icing conditions 9 apd 10, ^bie II. Results \bf calcu- 
lations of the rates of impin^ment aM eyaporat ion based the 
data are given in table IT. Ppr icing epnditim 9 tl» rate of water 
impingement was I.60 pounds per hour per foot span. rate of ' 
evaporation from the heated area was 0.44 ppuiid per hour per fppt 
span, leaving a calculated rate pf formation of ranbar.k of 1 . 1 6 p mmfl n 
per hour per foot of span. During icing condition 10 the calcinated" 
rate of iz^ingement was 1.79 poimds per 4010* per fppt and tim> rate of 
evaporation was 0.5I pound per hour per foot, resulting in a rate of 
formation of runback of 1.28 pounds per hour per fopt pf span, 
calculated rates of formation of runbapk (I.I6 and i,2& lb per hr, 
agree remarkably irell with the actual rate of formation (1,2 lb per 
hr, ft), .illiistrating the reliability of the ^pcedijre for caic^ 
latlng rate of evaporation. 

A shwt time priOT* to this test, the airfoil was subjected to 
a much less severe icing condition (condition 8, table It), during 
which eJJ. of the water intercepted was calculated to hayp been 
evaporated (icing condition 8, table IV). Photographs pf the test 
section verified the fact that no runback had formed. 

The foregoing ai^yses were based on the aesumptio4 that rempyal 
of the water striking the airfoil surface is effected by evaporation 
only, and that none of the water is dispersed by caI ~ metar^a j ' 

This is poi^istent with the results reported in reference IQ and 
25, It is believed that "blo^^ff" of water, as suggested in 
references 5 aM 6, does not occur. Also, it is believed there was 
no "boui^e-pff ” of the water drops striking the alrfpll stirfaces, as 
proposed in refereiwe 1^, At epeeds hi^er than thpse enccm[pass,ed 
by the scope of this Inyestl^gatl^, it is conceivable that mech^ical 
rempyal of the water by bounce:^ff could occiir, However, in vley pf 
the lack pf information pn this p^nonmnon, and since neglecting the 
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possltillty that water may he removed hy mechanical means tends to he 
more conservative in the thermal design, it is s\iggested that bounce— 
off he necp.ected in the design of wing thermal ice^prevention equip- 
3nent. 


Calculation of Heat Beqiiiremente for 
an HACA 0012 Airfoil 

Since it has heen demonstrated that the rate of evaporation of 
water from a heated wing can he calculated with reasonahle certainty, 
the rate of heat flow required to produce a particular rate of evapora- 
tion can he determined with equfO. dependability, provided the coeffi- 
cients of convective heat transfer are known. Using the equations 
and method piasented for calculating the rate of evaporation of water 
from a heated airfoil surface, a calculation was made to establish 

tho^extent-of^-heated area required_for-ice_pre_vention^Jn^pecif 

conditions of icing for the MCA 0012 airfoil, aesuming a particular 
heat— flow distribution. The conditions of calculation are as follows: 


Chord length 8 ft 

Pressure altitude .12,000 ft 

True airspeed 1?0 

Free-air teu^ierature .................... 20° F 

Llq\ild-water concentration .............. 0.3 gm/m° 

Mean-effective drop diameter 25 microns 


Drop-size distribution ......... ............F 

The procedure er 5 >loyed was to assume a reasonahle intensity and dis- 
tribution of total heat flow and then calculate the extent of heated 
area req^ziI*ed to evaporate all of the intercepted water. The method 
of solution will he outlined briefly in the followli^g paragraph. A 
detailed stejn-hy-step consideration of the problem showing all conq)U- 
tations is given in the appendix. 

First, the area, rate, and distribution of water lU 5 )lngement on 
the airfoil were calculated for the assumed conditions. Using the 
assumed distribution of total heat flow, the heat loss due to convec- 
tion for the particular conditions was then calculated. Since in an 
icing cloxid the presence of water on an airfoil surface causes pre- 
mature transition, for these calculations, transition was assumed to 
stsu^ at 5 percent s/c, 6ind the estimated form of turbulent heat- 
transfer coefficient shown in figure 25 was used. Calculations were 
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made for a number of chordwlse stations, and the results are given 
in figure 26, which shows the assumed heat-flow distribution and the 
calculated convective heat loss for one side of the airfoil. The 
rate of evaporation is represented by the area between the curves of 
convective heat loss and total heat flow, except in the region of 
water ii5)ingement, where the rate of evaporation is denoted by the 
area between the curve of heat loss duo to warming the intercepted 
water and the curve of heat loss due to convection, !Ihese areas 
act uall y give the rate of heat loss due to evaporation; however, by 
dividing the area value by Lg, the latent heat of vaporization, the 
rate of evaporation is obtained. The procedure, then, was to extend 
the total and convective cunres until the total rate of evaporation 
eqixaled the rate of water impingement. Extension of the heated area 
to l8 percent s/c was found to be adequ8''4d to ensure evaporation of 
all of the water intercepted. 

Several other calculations were made for the 0012 airfoil to 
determine the effects of altitude, air ten^jerature, and location of 
transition on -Uie requirements of heat flow and extent of heated 
area necessary to evaporate all the Intercepted water. The results 
of each of these calculations were conpared with the results of the 
calculations for the conditions previously specified. For each of 
the calculations, thei same total heat— flow distribution was assumed 
and the extent of heated area inquired to evaporate all the inter>- 
cepted water was calculated for each condition. 

To determine the effect of altitude on the heat reqtdzament, a 
comparative calculation was med.e for sea-^evel conditions with gi ?, 
other flight conditions as previously specified and with tto area 
and rate of water in^lngement the sam9 as at the 12,000-foot condi- 
tion, Kie results of this calculation are shown in figure 27, ^ich 
compares the relative convective heat losses at sea level 12,000 
feet. For the conditions at 12,000 feet extension of the heated 
area 'to l8 percent s/c was shown previously to be adequate to 
ensure evaporation of all of the water Intercepted, At sea level it 
woxxld be necessary to extend the heated area to 26 percent s/c for 
evaporation of all the water intercepted. The curves of figure 27 
can also be used to determliie the amount of Increase necessary in the 
total heat flow if aJ.1 the water Is to be evaporated in an area 
forward of a specified chord point. For exas^le, assume that the 
extent of heated region for the curves of figure 27 is limited to 
l8 percent s/c. At 12,000 feet all of the water woxild be evapo- 
rated, as has been previously mentioned. At sea level, however, 
some of the water would not have been evaporated. By measurement of 
the areas of figure 27 it can be shown that the total heat flow 
required to evaporate all the water within the area from 0 to l8 
percent s/c at sea level is approximately 10 percent greater than 
the amount required at 12,000 feet. The increase in heat requirement 
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with decrease In altitude is due to the fact that the rate of evapo- 
ration of water decreases as altitude is decreased, "because of the 
decrease in the evaporative factor X (equation (23) )• Since the 
convective heat— transfer coefficient increases with decrease in alti- 
tude, due to the increase in air density, it mi^t he expected that 
the rate of evaporation would he increased with decrease in altitude, 
because the rate of evaporation is directly proportional to the con- 
vective coefficient (equation (21)). However, the increase in the 
rate -of evaporation is more than conq)ensated hy the increase in con- 
vective heat loss, and the rate of evaporation, for a fixed total 
heat flow, actually hecomes less with decrease in altitude. Appar- 
ently, then, ' airfoil thermal ice— prevention equipment in which the 
heat flow is fixed, such as electrical systems, should he designed 
for the minimum altitude at which the airplane is expected to 
encounter icing. However, if the airplane is designed to utilize 
some form of air-heated system, the performance of which prohahly 

-will- decrease- Afith -iiicrease_in_alh.ltude,_ -toe ^laximum altitu de at 

which icing is expected to he encountered should also he investi- 
gated. 


To determine the effect of air tenqierature on the heat require- 
ment, a calculation was made of the convective heat loss at 0° f 
free-air temperature and is conq>ared in figure 27 with the convective 
heat loss at 20° F. In the calculation with the free-air temperature 
at 0° F, it was determined that the surface temperature dropped to ■ 
freezing at 24 percent s/c before 'all the water on the surface was 
evaporated. However, the total heat flow required to evaporate all 
the water within the area from 0 to I8 percent s/c with the air 
temperature at 0° F is approximately only 15 percent greater than the 
amount required at 20° F. Althou^ this is an appreciable increase 
in the heat requirement, it is considerably less than that necessary 
for a similar change in conditions for ice-prevention equipnent 
designed on the basis of maintaining the surface temperature ' Just 
above freezing# such as for the case of windshields. (See infer- 
ence 21.) It appears, then, that a wing thermal system which has 
been designed for a relatively hi^ air temperatxxre will be capable 
of. ice prevention at low air tenq^ratures in icing conditions nearly 
as severe as those upon which the design was based. Of course, the 
system is more subject to failure through the possibility of the 
surface ten^^rature falling below freezing in the low air-tenq)erature 
conditions, but in general, the surface temperatxires required for 
evaporation of all iDq)inging water in the relatively small heated 
area of the leading edge will be sufficiently high to obviate this 
possibility.. 

To establish the effect of the location of transition on the heat 
requirement, a calculation was made of the convective heat loss, 
assuming laminar flow exists throu^piout tho heated area. For this 
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ealculatloa> the SBaeured 76dues of cotivectlTe heat-transfer coeffl— " 
elent shoim in figure 2^ were used. The conrective heat loss for 
laminar flov is conqpared in figure 27 with the convective heat loss^ 
assuming transition started at ^ percent s/c. In the case of complete 
laminar flov, it vould he necessary to heat only to lU percent s/c 
to obtain evaporation of all the water. The total heat flow required 
to evaporate all the water within the area from 0 to Ih percent s/c 
with transition at 5 peircont s/c is approximately 10 percent greater 
■Uian the amorunt of heat requli^d if laminar flow prevails. Apparently, 
the location of transition moves forward in conditions of icing, even 
in the presence of a favorable pressure gradient, to a point idiere a 
strong favorable iressure gradient is encountered (figs. 11 and l6). 

As was stated previously, the location of transition is believed to 
fluctuate, probably over a considerable distance. It is stiggested 
that forward movement of transition to a point close to the leading 
edge of the wing be assumed in the design of thermal ice— prevention 
equlpnsnt, especially in view of the fact that a greater amount of 
heat is required for the turbulent-flow condition. 

Prom a conqirehenslve study Of the results shown in figure 27, 
some general conclusions can be reached. It is apparent that eift of 
the area of droplet Impingement, the efficiency of removal of water 
by evaporation decreases rapidly. The reason for the decrease in 
efficiency is that (xily partial wetness prevails aft of the area of 
isplngement, while the area of Ingpingement is entirely wet. This 
indicates tMt the larger the portion of the total amount of water 
Intercepted that is evaporated in the area of interception, the 
greater the efficiency of the. thermal system becomss. The rate of 
evaporation of water is the determining factor in the efficiency of 
a wing thermal ice— prevention system. Only the heat that is dissi- 
pated in evaporation is used to advantage, Ihe heat lost by conveo- 
tlon only warns the air. Thus, the conclusion is drawn that the 
heating should be concentrated as much as possible in the leading 
edge of a wing, in the area of drop lsq>ingement, if an efficient 
thermal system is to be obtained. 


Calculations for the C-46 Wing Thermal System 
in Maximum Continuous Icing Conditions 

V 

An analysis of the 0-46 airplane wing thermal ice— prevention 
system for the upper sinrface at wing station 157 was made in an 
effort to determine vbethsr the thermal system could cope with the 
mailmum continuous icing conditions given previously at the beginning 
of this discussion. The assumed icing and flight conditions for the 
calculations are as follows: 
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Condition A 

Condition B 

Altitude (ft) 

6000 

6000 

True airspeed (n^h) 

180 

180 

Liquid water 
content (gm/m®) 

0.8 

0.5 

Mean-effective 
drop size (microns) 

15 

25 

Eree-alr 
ten^ierature (°P) 

20 

20 


The hea-b-flow dietribution at station 157 was estimated, based on 
data presented In references 4 and 26, and Is shown in figure 28. 

The rates of water ia^ilngement ter the two Icing conditions assuming 
a ”C" type drop-size distribution were calculated for the leading- 
edge cylinder of the airfoil section using the data presented in 
reference 17. Curves of distribution of water impingement, for the 
upper surface, for the two cases are given in figure 29, These were 
constructed using the data from reference I7 for each drop size in 
the distributions. The value of Mg, obtained from equation (9), 
for Condition A is O.65 pound per ho\ir per foot span, while the vedue 
of Kfe for Condition B is 1.39 pounds per hour per foot span. As in 
the calculations presented ia?eviously, the rate of evaporation of 
water from the surface was determined by calculating the heat loss 
due to convection for the two conditions. Values of convective heat- 
transfer coefficient were taken from figure 30, which shows the values 
msasmed during fli^t in clear air using an electrically heated glove 
And the estimated convective coefficients for icing conditions, \dien 
transition moves forward. Uie estimated values were used in the cal- 
culations. Ihe congouted curve of convective heat loss for Condition 
A is shown in figure 28. Besults of the calcxilations of rate of 
evaporation for the two conditions Indicated that sufficient heat 
' was supplied to the upper wing siirface to evaporate all of the water 
intercepted. For Condition A the wing was shown to be capable of 
evaporating O.9O pound per hour per foot span, indicating that the 
liquid— water concentration .could attain a value of at least 1.1 grams 
per cubic meter at a mean-effective drop size of 15 microns before 
runback would form. The rate of evaporation for Condition B was 
calculated to be I.35 pounds per hour per foot span, sviggestlng that 
the liquid— water concentration of O.5 gram per cubic meter is the 
limiting condition at 25 microns. 
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A fxirther analysis was made of the upper wing surface at station 
157 using the values of Conditicai A, excepting that a free-air tempera- 
tvire of 0° F was assumed* The curve of convective heat loss for this 
condition is shown in figure 28* Ifoder this condition, calciilations 
revealed the wing would he able to evaporate 0*70 pound per hour per 
foot span. The rate of iD5)ingeinent, as before, was 0,65 pound per 
hour per foot span. 

These calculations substantiate the general observations of the 
successful operation of the C-h 6 wing themal ice— prevention system. 

The absence of runback on the wing upper surfaces dtiring a great many 
of the icing flights indicates the adeq.uacy of the thermal design. 


Selection of Conditions for Design 

In selecting values of drop size, liquid-water concentration, 
air temperature, and altitude for the design of thermal ice— prevention 
equlpnent, a combination of these variables normally occxirrlng in 
nature should be chosen such as to require the hipest rate of heating. 
As stated previously, conditions of maxlmtmi continuous icing are 
believed to form a good basis for design. It is of Interest to inves- 
tigate the effect of different possible combinations of the variables 
of drop size, liquid— water content, air tenq)erature, and altitude on 
the heat requirement for ice prevention for the maximum continuous 
conditions given in the table in the first part of this discussion. 

The effect of an increase in the size of drops in an Icing condi- 
tion is to increase the collection efficiency of the airfoil, thereby 
increasing both the rate at which water is Intercepted and the area 
of inqjin^ement. An increase in the liquid— water concentration of 
the air causes a proportional Increase in the amount of water inter- 
cepted, for a given drop size. Since all of the water striking the 
wing must be evaporated to avoid the formation of ice, an Increase 
in the rate of water Interception will cause an Increase in the heat 
requirement. Fortunately, a relation between water concentration 
and drop size appears to exist in icing clouds, and the existence of 
very large drops generally is accompanied by a small concentration 
of liquid water (reference I3). The selection of a combination of 
drop size aind water concentration should be such as to produce the 
highest rate of inq)ingement. It was shown previously that eui Increase 
in the drop size produces a greater Increase in the rate of water 
interception than a proportional increase in the water concentration. 
For this reason, the maximum continuous icing condition of 25 microns, 
mean-effective drop size, and 0.5 gram per cubic meter, liquid- water 
concentration, generally will result in a more rapid rate of water 
inq)lngement than the condition of I5 microns and 0.8 gram per cubic 
meter. 
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A decrease in free-air temperature, vhile Increasing the heat 
requirement for thermal ice prevention, is acconqianied by a decrease 
in the liquid-isater concentration (reference lU), idiich causes a 
proportional decrease in the rate of water in 5 )ingeinent, for the same 
drop size. The sizes of drops existing at low air ten 5 )eratxjree (O® P) 
In icing conditions tend to be only sli^tly smaller than those at 
hi^er air teinperatures (reference 13 )j therefore, a selection of 
air tonq)erattjre for design will be determined by the combination of 
air tenq)erature and water concentration (and corresponding drop size) 
to produce the highest heat requirement. It will be found, genereilly, 
that the rate of heating required to evaporate the larger quantities 
of water at the higher air temperatures is greater than the heat 
needed for ice prevention at the lower teiiperatures. However, low 
aij>-tenperature conditions should be investigated to ascertain that 
the tempexature of the heated surface will not fall below freezing. 

TSeie^appeara to-be-no-relatlon-between^altitud e and the dro p 
size or liqtiid-water concentration of icing conditions. (See refeiv^ 
ence I 3 .) Therefore, the altitude at which the heat requirement is 
greatest should be chosen. The TniniTmim altitude of operation was 
shown previously to produce the hlg^st heat acquirement for wing 
thexmal ice pacvention. However, as was foavnerly suggested, if the 
alirplaaie is designed to utilize some fom of aia>-beated system, the 
mftTiimnn altitude at which icing is eacpected to be encounteacd should 
also be investigated. 


C0HCLH5ICHS 

Prom the foicgoing discussion, it is concluded that the eartent 
of knowledge on the meteorology of icing, the impingement of water 
drops on alai^foil surfaces, and the processes of heat tiransfer and 
evaporation from a wetted aiarfoil suarface has been Ixacacased to a 
point wheac the design of heated wings on a fundameiatal, wet-air 
basis now can be uiadeartaken with acasonable ceartainty. In addition 
to this genead conclusion, the following coiacluslons are drawn, 
based on test data and aaaalytlcal studies of the processes of heat 
tacnsfer evaporation from a heated wing: 

1. The heat should be concentacted as much as possible in the 
leading-edge region of the wing in the area of wateavdrop impinge- 
anent, if an efficient thearmal system is to be obtaiiaed. 

2. An iaacacase in altitude, for the same acte aaad aarea of water 
impingement on a wing for the same coxaditlons of tarue airspeed 
and face-air temperatuac, decacases the heat requlacment for theavnal 
ice preventlm. 
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3« A vlng thermal ice— preyentlon system which has heen designed 
to evaporate all impinging water in the leading-edge region for a 
relatively hi^ free-air ten5>era.txire (20° F) will be capable of ice 
prevention at low air tenqwratures (0° F) in icing conditions nearly 
as severe as those upon which the design was based. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX 

Calculation of Extent of Heated Area'Eequlred 
for NACA 0012 Airfoil 

The detailed, step-by-step calculations for establishing the 
extent of heated area required for ice prevention on an NACA 0012 
airfoil in specified conditions of icing are presented in this appen- 
dix. It is believed that the general procedure outlined herein will 
be applicable for the design of most wing thermal ice— prevention 
equipment. 

The calculations wore made for one side only of the airfoil. 

The assumed flight and meteorological conditions used in the calcu- 


lations are as follows: 

Pressure altitude ......... ... 12,000 ft 

True airspeed ITO s^ih 

Free-elr temperatiure. 20° F 

Liquid— water concentx^tlon. 0.3 gm/m° 

Mean-effective drop diameter. 23 microns 

Drop-size distribution. ..................... E 


The chord length of the airfoil was taken as 8 feet. 

Step 1 .— Calculate area, rate, and distribution of water inter- 
ception. The area of interception is determined by the largest 
droplets present in the cloud. For the case of an "E" type drop- 
size distribution (reference I7), the largest drops will be 2. 71 
times larger than the mean-effective size, or 
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Maximum drop diameter = 2,71 x 25 = 68 microns 

The K— value corresponding to a drop diameter of 68 microns for the 
assiamed flight conditions vas cedculated vising equation (5) , Thus, 

K = a ( l4xii28>ao[l Y ( 250x8x0,05^ , 0.465 

9 ^0.053'' N 8 / \ i.l6xlor® ^ 


From figure 3> "the efficiency of ioq)ingement E for this K—veilue is 
54 percent. Using equation (10), the starting ordinate of the 68 - 
micron-diameter drop which Just hits the surface is 


y°limit " ® 


or 


y°limlt _ ymax 

c c 


Since the airfoil is 12 percent thick. 


, 0.06 

c 


and 

^°limit = 0.54 x 0.06 = O.O 32 
c 

Using the broken curve in figvu^ 1, the area of inq)ingeinent was found 
to be 10.8 percent s/c. 

'Hie distribution of water-impingement rate over the airfoil 
surface was calculated using eqxiatlon ( 8 ), The IndlviduEil rates of 
Inqiingement for each of the drop sizes in the assvoned distribution 
were calculated for various points along the surface. This was 
accomplished by computing the value of K for each of the drop sizes 


k2 
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in the distribution using equation (5), determining the values of C, 
for the corresponding K-value, at various points along the svirface, 
then evaluating the expression VnC. By adding these individual 
iDq)ingement rates for each point on the surface, the resulting dis- 
tribution of vater-in5>ingement rate over the surface was obtained. 

The above calculations were made using a tabular form of conq)utation, 
as Illustrated in table V. Figure I7 shows the resulting distribu- 
tion of impingement. The total rate of water impingement Mg was 
calculated, using equation (9) , by measuring the area under the cui’ve 
shovm in figure I7. The value of Mg was ceaciilated to be 2.1 
pounds per hour per foot span. 

Step 2 .— Establish the distribution of heat flow from the 
surface. The distribution of heat flow will depend on the type of 
ice— prevention system to be used. If an electrical system is 
planned, the distribution and Intensity, once set, will remain 
xmchanged regardless of variations in flight and meteorological 
conditions. On the other hand, if the system is to be designed to 
utilize heated air or some other fluid, the distribution of heat flow 
from the surface will depend upon the characteristics of the internal 
flow of the fluid as well as the conditions affecting the external 
heat transfer. If such a system is to be used, calculation of the 
heat-flow distribution will be rather con5>lex, and it is believed 
that assuming a distribution will provide a good basis for starting 
the calculations for design. 

The heat— flow distribution and intensity used in these calcu- 
lations was estimated, based on dfta presented in references 4 and 
26, to be the heat— flow distribution and intensity of the thermal 
ice— prevention system for the wing of the C—kS airplane (refer- 
ence 20) at the 8-foot chord station. This distribution, which Is 
believed to be representative of a probable thermal system, is shown 
in figure 26, 

Step 3 »— Determine the values of convective heat— transfer 

coefficient. The values of measured convective heat— transfer coeffi- 
cient with the estimated form of the turbulent coefficients shown in 
figure 25 were used. 

Step 4 .— Calculate values of surface temperature in the area of 
water impingement, using equation ( 2 k) ^ such that the values of q at 
any point are equal to the assumed heat flow. (See fig. 26.) 

Step From the values of surface temperature calculated in 

step k, compute the convective hea't loss using equation (12). The 
curve of convective heat loss is plotted in figure 26. 
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Step 6 .— Calculate the rate of flow of water aft of the region 
of water interception. This was done hy measuring the area in the 
region of water impingement between the convective heat— flow curve and 
the curve denoting heat flow to inq)lnglng water (fig. 26) to obtain 
the rate of heat dissipation due to evaporation. The rate of evapora- 
tion was coD5)uted from equation (20), and was subtracted from the rate 
of impingement to give the rate of water flow aft of the region of water 
interception. The rate of evaporation was calcxilated to be 1.6 pounds 
per hour per foot span. Subti^cting this value from the rate of water 
striking the surface, 2.1 pounds per hour, foot, the rate of flow of 
water aft, then, is 0.5 pound per hour per foot span. 

Step 7 »— Determine the wetness fraction and make the proper 
modification to the evaporative factor. Using the curve shown in 
figure 21, the wetness fraction for a water flow rate of 0.5 pound 
per hour, foot is 30 percent. It is suggested that the values of 

de^^^of " iretMs's~ glven~in-f igure— 21-be--used-only--to_the _nearest. 

10 percent, since more precise usage is considered to be unwarranted. 

The evaporative factor X was then modified by the 30-P©rcent 
wetness fraction, so that equation (23) becomes 



For these calculations, the value of Pi was taken as the free— stream 
static pressure, so that 


X = 1 + 1.77 


tg-t 


2c 

Ok 


) 


Step 8 .— Calculate values of surface ten5)erature aft of the area 
of water iBq)ingement xising equation (25) and the revised value of X 
such that the values of q at any point are equal to the assumed heat 
flow. (See fig. 26.) 

Step 9 .— From the values of surface tenqjerature calcxilated in 
step 8, compute the convective heat loss, using equation (12). The 
curve of convective heat loss is shown in figure 26. 

Step 10 .— Extend the calculation of convective heat loss xuitll 
the total rate of evaporation, denoted by the area between the cin*ves 
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of convection and total heat flow (except in the region of in 5 )ingenient), 
equals the total rate of water i 25 >inge 3 nent. The rate of evaporation 
was con 5 )uted from the area hetween the two curves using equation (20), 
For the case of the 0012 airfoil, the extent of heated region required 
for evaporation of 2.1 pounds per hour, foot span was calcxxlated to he 
to l8 percent s/c, which is equivalent to l6,5 percent chord. 

It shoxild he noted that the extent of the heated region can he 
decreased hy increasing the intensity of the total heat-flow distri- 
bution arifl re-calculating the required extent of heated area. 
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KACA TN No. 1 ^ 4-72 


TAEIE m.- CCMPARISQR OF CALCX1LA!IED BATES OF WAOEB 

iMPUfcacMssT aud evAfqratiqn over the leading edge 

OF THE NACA 0012 EIEdBICALLT SEATED AIBFOIL 
MODEL FOR ICING CONDITIONS OF TABLE I 





Calculated 

Cedculated 




rate of 

rate of 

Icing 

condi- 

tion 

ni^t 

number 

Pacific standard 
time 

vater 
ii^inge- 
ment, % 
[ lb/(hr) 
(ft. span)] 

eater 
evapora- 
tion, Ws 
[ lb/(hr) 
(ft. span)] 

1 

39 

2:03 to 2:05 

0.30 

0.17 

2 

39 

2:13 to 2:16 

.35 

.in 

3 

39 

2:19 to 2:22 

.27 

.36 

k 

39 

2:23 to 2:26 

.03 

.04 

5 

39 

2:28 to 2:31 

.19 

.19 

6 

1^3 

12:27 to 12:29 

.61 

.62 

7 

^^9 

1:25 to 1:26 

i.in 

.24* 


^K)Dly leadlne-edgB region heated. Bvinback fonaed. / 

.... ,/ 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS. 
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TABLE 17.- C(M*ABISQH OF CALCULATED HA3ES OF WA3EB 
IHPIHGS!ME37r AHD E7APC!RA!ITQn OVER THE LEADIHU 
EDGS OF THE HACA 6^,2-016 ELECTRICALLY 
TTT^TCD AIRFOIL HCEDEL FOR ICIHG 
CONDITIONS OF TABLE II 


Iciixg 

_COT|U.— 

tion 

Flight 

-nimiber- 

Pacific Standard 
time 

Calculated 
rate of 
watex'-drop 
impinge— 

- ment, — % 
[ l‘b/(hr) 
(ft.span)] 

Calculated 
rate of 
water 
evapora- 
__tlon, Wg 
[ lh/(hr) 
(ft.span)] 

1 

100 

3:15 to 3:19 

0.47 

0,46 

2 

102 

1:30 to 1:35 

.94 

.96 

3 

105 

11:01 to H:06 

.52 

.55 

h 

105 

11:1H to 11:45 

.56 

.54 

5 

105 

2:32 to 2:36 

.97 

1.05 

6 

105 

3:07 to 3:12 

1.01 

1.06 

T 

105 

3:18 -to 3:23 

.33 

.38 

8 

111 

12:32 to 12:36 

.18 

0 

OJ 

• 

9 

111 

1:08 to 1:12 

1.60 

.44* 

10 

111 

1:15 to 1:19 

1.79 

.51* 


♦Only leading-edge region heated, Eunhack formed. 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS. 
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><''■ o 



« 


Figure 5*— C-46 test airplane as flown during the winter of 19^6-47 showing the position 
airfoil models mounted on the fuselage. 
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(a) NACA 0012 section moimted on C— ^ fuselage 
for the 1945-*^ fll^t tests. 



-CJSs 


NACA 


A-11516 


(h) NACA 65,2-016 section mounted on C— 46 fuselage 
for 1946-47 fll^t tests. 

Figure 6.— Electrically heated airfoil models used to obtain data 
in natural— icing conditions. 
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8 


7 ^ 


4. 


/. Surface thermocou0es' ^ 

2 Aluminum 'skin 

3 Plastic— impregnated fabric ; 

4. Electrical resistance heating strips 

5. Plastic- impregnated fabric 

6. Plastic base 

7 AirfoH structure 

8 Leads from surface thermocouples '' 

9 Electrically heated guard sections 

to Electricaf leads from resistance heafmg 
//. Viewing blister on fuselage 
12. Camera 


Figure 7.- Cut-eway view of the RACA 65,2-016 electrlceilly 
heated airfoil model showing construction details. 
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Figure 9 .— l^pical rnintack fonnation obtained on the HACA 0012 airfoil 
model with only leading-edge region heated. 
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Figure 13.— Runback formation obtained on the electrically heated 
NACA 65,2-016 airfoil model during Icing conditions 9 10* 

table II. 
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Figure 15*— Typical data records showing area of drop Impingement on the KACA 65,2-Ol6 airfoil 
model, and traces of water flow aft of the region of Impingement. 
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